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Extremely Red Galaxies: age and dust degeneracy solved ? 
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ABSTRACT 

The extremely red galaxies (ERGs) are defined in terms of their very red optical-to- 
near IR colours (as R-K>5 or I-K>4). Originally this selection was aimed at selecting 
old (> 1 Gyr) passively evolving elliptical galaxies at intermediate redshift (l<z<2), 
but it was soon discovered that young star-forming dusty galaxies can show similar 
colours and therefore be selected in the same surveys. It is crucial to distinguish 
between these two populations because they have very different consequences on the 
models of galaxy formation. Here we show that old ellipticals and dusty starburst are 
expected to show different colours in the (I-K) vs. (J-K) diagram for redshift range 
l<z<2, providing thus a useful tool to classify ERGs in large samples up to K<20. 
This is mainly due to the fact that old galaxies at these redshifts have a strong 4000A 
break at A < 1.2/im (J band), while dusty galaxies show a smoother spectral energy 
distributions and therefore redder J-K colours. We discuss this difference in detail 
both in the framework of the stellar population synthesis models and by using observed 
spectra. The selection criterion is also compared with the properties of ERGs of known 
nature. We also show that this colour selection criterion is also useful to separate the 
ERGs from brown dwarf stars showing similar optical-to-IR colours. 

Key words: galaxies: elliptical and lenticular, CD - galaxies: formation - galax- 
ies: fundamental parameters - galaxies: starburst - galaxies: evolution - cosmology: 
miscellaneous. 



1 INTRODUCTION 

The extremely red galaxies (ERGs) discovered in deep near- 
infrared (IR) and optical surveys (e.g., McCarthy et al. 
1992, Hu & Ridgway 1994) are among the most intriguing 
objects in the modern cosmology and their nature is still 
controversial. They are defined in terms of their very red 
optical/near-IR colours (R — K > 5 or I — K > 4), are very 
rare at bright K magnitudes while their density approaches 
0.5 ±0.1 arcmin" 2 at K=20 (McCracken et al. 2000). Such 
very red colours can be explained by two mainly oppos- 
ing scenarios: 1) ERGs can be starburst galaxies hidden by 
large amounts of dust, or 2) can be high redshift (z > 1) 
old ellipticals with intrinsically red spectral energy distri- 
butions (SEDs) and large positive k-corrections. To know 
which population is dominant between the ERGs is a cru- 
cial test for the models of galaxy formation: the presence of a 
wide-spread population of high redshift old ellipticals would 
point towards a very early formation of massive galaxies 
in a "monolithic" scenario (Eggen, Lynden-Bell & Sandage 
1962, Larson 1975), while the presence of numerous dusty 
starbursts at intermediate redshift is better fitted by the 
hierarchical scenarios of galaxy formation (White & Frenk 



1991, Kauffmann, Chariot & White 1996). Both these pop- 
ulations would have a deep impact on the actual knowledge 
of the cosmic history of star formation. 

The two populations can in principle be distinguished 
by several observations: 1) high signal-to-noise spectra can 
show the presence of emission lines, revealing ongoing star 
formation activity, or, on the contrary, sharp spectral breaks 
due to old stellar populations; 2) sub-mm observations can 
detect the far-IR emission from the hot dust heated by the 
starburst; 3) high resolution optical and near-IR images can 
reveal the morphology, which is expected to be regular and 
centrally peaked for old ellipticals and irregular, disturbed 
or even double for dusty starbursts. In practice these ob- 
servations help only in a limited number of cases: 1) the 
faintness of these objects in the near-IR and especially in 
the optical (typical magnitudes are K=19 and R=25) makes 
the spectroscopic observations very difficult, and many Keck 
and VLT nights have produced redshifts for just few ob- 
jects (see below); 2) the best current sub-mm instrument, 
SCUBA (Holland et al. 1999), can still detect the far-IR 
flux only from exceptionally bright objects (e.g., Cimatti et 
al. 1999) , while the vast majority of the population remains 
below the detection threshold; 3) HST (WFPC2 /NICMOS) 
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Figure 1. Representative spectra from an old elliptical (thick 
line) and a dusty starburst (thin line) at z=1.5. The position of 
the I, J and K filters are shown. Both spectra have been computed 
using the Bruzual & Chariot (1999) GISSEL model. The thick line 
is a simple stellar population (SSP) with an age of 15 Gyr, while 
the thin line shows the spectrum resulting from a constant star 
formation rate over 1 Gyr when reddened to E(B-V)=0.8 by a 
screen of dust described by the SMC law. It is apparent as, while 
the I-K colour selects both class of objects, the J-K colour can 
distinguish between them. 



images have been obtained only for a limited number of 
objects, and also in these cases a detail optical morpholog- 
ical study is hampered by the faintness of these galaxies. 
Therefore it would be useful to discover a classification tool 
applicable to large samples of ERGs. 

A few ERGs were studied in great detail and can be 
used as benchmarks for a classification. As expected, both 
dusty and old galaxies contribute to the ERGs population: 
HR10 (Hu & Ridgway 1994), a colour-selected galaxy at 
Z—1A4 (Graham & Dey 1996), shows a disturbed and elon- 
gated WFPC2 morphology (Dey et al. 1999) and was also 
detected by SCUBA in the sub-mm (Cimatti et al. 1998, 
Dey et al. 1999), making its classification as a dusty star- 
burst very robust; 53W091 (Spinrad et al. 1997), a faint 
radio galaxy at z=1.55, whose spectrum shows a break due 
to a dominant stellar population of a few Gyr old (Heap et 
al. 1998); 6 galaxies in the core of the X-ray selected clus- 
ter RXJ0848.9+4452 (Rosati et al. 1999) at redshift 1.26 
show absorption features typical of evolved ellipticals and 
3 ERGs with K < 20 (Cohen et al. 1998) are found to be 
absorption line galaxies at redshift 1 < z < 1.23. Some 
ERGs still without redshift show no evidence for strong 
emission lines (Cimatti et al. 1999, Cohen et al. 1999), two of 
them have however SEDs that require strong dust reddening. 
Other ERGs can be classified by their HST optical morphol- 
ogy (Pozzetti et al. in preparation) and surface brightness 
profiles (Moriondo et al. 2000). Finally some faint sub-mm 
sources have been identified with ERGs (Smail et al. 1999). 

All the ERGs with spectroscopic redshifts fall in the 
range 1 < z < 2. Moreover, most of the other ERGs dis- 
covered by the current surveys are expected to fall in this 
redshift range because of their optical and near-IR colours 
and magnitudes. Indeed L* objects at higher redshifts are 
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Figure 2. Colors predicted for ellipticals (thick line) and star- 
burst (thin line). The dotted lines represent the colours for z<l 
and z>2, while the solid lines show the redshift range Kz<2. 
The dashed lines indicate the ERGs selection criterion (I-K>4) 
and separate the region where the two populations of ERGs are 
expected to fall. Stars and dots show the stars in the samples 
by Pickles (1998) and the brown dwarfs by Legget et al. (1998), 
respectively. The I-K colours of these stars are similar to those 
of the ERGs, but the two populations show different J-K average 
colours, the brown dwarfs being bluer. 

expected to fall below the magnitude threshold of the cur- 
rent surveys (K < 20), while objects at lower redshift tent 
to drop out of the colour selection for ERGs (R — K > 5 or 
/ — K > 4) because of smaller k-corrections. For these rea- 
sons we have studied a photometric way to distinguish be- 
tween the two populations in the redshift range l<z<2. This 
recipe is based on optical/near-IR broad-band data only to 
be easily applied to large samples of ERGs; it is derived 
by using numerous stellar synthesis models and observed 
spectra. The prediction can be compared to the properties 
of observed ERGs for which we can assign a classification 
based on morphological, spectral or sub-mm information. 



2 THE I-K VS. J-K COLOUR PLANE 

While old stellar populations are characterized by sharp 
breaks, especially around 4000 A, young dusty galaxies show 
shallower SEDs in the optical range because dust extinction 
does not produce sharp breaks. This is shown in Fig. 1 where 
the representative spectra of the two classes of objects are 
shown at z — 1.5. This feature allows us to clearly sepa- 
rate the two classes of objects using a I-K vs. J-K diagram, 
where the former colour measures how red an object is, while 
the latter measures the steepness of its SED. Daddi (pri- 
vate communication) propose a similar classification scheme 
based on the J-K colour only. 

In Fig. 2 we show the I-K vs. J-K plane when local 
galaxies are used as templates of high redshift objects: for 
the ellipticals we have used the average of 5 local ellipticals 
by Mannucci et al. (in preparation), while we have used M82 
(see Silva et al. 1998) as a typical dusty starburst. When 
redshifted between z = 1 and z = 2, these observed spectra 
fall in different regions of the plane and for a given I-K they 
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are separated by about 0.3 in J-K. In this diagram, as in 
the following ones, we use a Cousins' I filter and Bessel & 
Brett (1988) J- and K- filters. In this figure we have also 
plotted some cold stars and brown dwarfs which are known 
to contaminate the ERGs samples selected by the optical- 
to-IR colour criterion, especially at faint flux levels when the 
morphological classification becomes weaker. When the J-K 
colours is considered, the two populations seem to clearly 
separate (see also Fig. 4), at least when the cold stars in the 
Pickles (1998) library and the brown dwarf in the Legget et 
al. (1998) are considered. Similar results are also obtained 
by using models and observations by D'Antona et al. (2000) 
and Burrows et al. (1997). The use of the J-K colour can 
therefore remove most of the foreground red stars from the 
ERGs samples. 

In Fig. 3 we show that this separation is predicted by 
a large range of models of old elliptical and dusty starburst. 
We have explored this range by studying how the results 
depend on the particular spectral synthesis model used, the 
metallicity of the stellar population, the assumed history 
of star formation, the Initial Mass Function (IMF), the as- 
sumed redshift of formation and the properties and the dis- 
tribution of dust. 

The various panels of Fig. 3 show the results for l<z<2. 
In all the cases the dashed line from Fig. 1 (see eq. 1) sepa- 
rates the two populations, ellipticals on the left and starbust 
on the right. We used a cosmology with Ho = 50 km s^ 1 
Mpc -1 and — 0.1, but this choice is not critical for the 
results. 

• Panel (a): The basic model. 

The ellipticals are represented by simple stellar population 
(SSP) of the Bruzual & Chariot (1999) evolutionary model 
with solar metallicity (Z=0.02), Salpeter IMF between 0.1 
and 125 Mq and three values of the redshift of formation, 
2/ = 3,4,6. The starbursts are reproduced with constant 
star formation rate, z/=2.5 and 6, and an extinction by dust 
with E(B — V)—0.5 and 0.8, assuming the Small Magellanic 
Cloud (SMC) extinction law (Prevot et al. 1984), and a uni- 
form screen external to the emitting region. These models 
reproduce the local ellipticals and starburst very well. De- 
spite the spectra used for Fig. 2, these spectra take into 
account the evolution of the galaxies with look-back time. 

• Panel (b): dust properties and distribution. 

We explored the effect of different dust extinction laws and 
distributions in starbursts, in addition to the SMC law used 
in the other panels. Firstly we have used the Calzetti (1997) 
law derived using starburst galaxies, assuming, as for the 
SMC one, that dust is external of the emitting region. This 
may be not appropriate to reproduce starburst galaxies 
where dust in supposed to be very inhomogeneously dis- 
tributed and mixed with the emitting gas. Therefore we also 
used the extinction laws by Silva et al. (1998), Ferrara et al. 
(1999) and Cimatti et al. (1997), based on a number of dif- 
ferent assumptions for dust spatial distribution and optical 
depth, including both absorption and scattering by dust. 

• Panel (c): metallicity effects 

The same models of panel (a) are plotted but with two dif- 
ferent metallicities, Z=0.004 (0.2 solar, bluer curves) and 
Z=0.05 (2.5 solar, redder curves). The predicted colours be- 
come bluer and redder respectively, but each class remains 
inside its region of the diagram. 



• Panel (d): IMF effects 

Here a Scalo IMF is used, which is steeper at high masses 
and shallower at low masses than the Salpeter IMF of panel 
(a). The resulting effect is a prediction of colours redder than 
in the basic model. 

• Panel (e): Star formation histories 

We also changed the history of star formation by assuming a 
small residual star formation for elliptical galaxies described 
by an exponentially declining SFR with a short e-fold time 
of 0.3 Gyr, and for dusty starburst an e-fold time of 0.1 
Gyrs, at an age of 0.1 Gyr and for E{B - V) = 0.5, 0.8. 

• Panel (f): other models 

Independently from cosmology we plot no-evolution mod- 
els for elliptical (SSP at age=l, 5 and 15 Gyrs) and dusty 
starburst (constant SFR and age=0.1,l and 15 Gyrs and 
E(B-V)=0.5,0.8). We also want to explore how the predic- 
tion change when different spectrophotometric models are 
used. We use Barbara & Poggianti (1997) models (dotted 
lines), SSP for elliptical galaxies at age of 1, 15 Gyrs and 
constant star formation models at age=l, 15 Gyrs with an 
extinction of E(B-V)=0.5,0.8 for starbursts. 

In conclusion, by changing the various parameters of the 
stellar synthesis models it is poSsible to obtain quite differ- 
ent spectra, but the separation between the two populations 
remains a valid prediction. Since some models fall near the 
separation line, some degree of contamination between the 
two class could be present near it. 



3 ERGS OF KNOWN NATURE 

These predictions can be compared with the observed prop- 
erties of ERGs with available I, J and K magnitudes and 
whose nature can be independently determined from spec- 
tra, morphology or the presence of a detectable sub-mm flux. 
In total we found 24 objects, 6 starbursts and 18 ellipticals, 
that can be placed on the I-K vs J-K diagram . The re- 
sults are shown in Fig. 4 (see the caption for references). 
Although the number of objects is still small we found that 
all the objects fall in the part of the diagram predicted by 
the models. Some objects of both classes fall near the sep- 
aration line: the classification in this range should be used 
with caution especially in a sample with poor photometry. 



4 THE SELECTION CRITERIA 

The dashed line in Figures 2, 3 and 4 which separate the two 
populations for l<z<2 in the I-K vs. J-K plane corresponds 
to: 

(J - K) = 0.36 • (7 - K) + 0.46 and / - K > 4 (1) 

and is consistent with the redshifted spectra of local objects 
(Fig. 2), the predictions of the spectral synthesis models 
(Fig. 3) and the ERGs of known nature (Fig. 4). 

The filter R (Cousins) could also be used instead of I to 
separate the different ERGs populations at similar redshift. 
By arguments very similar to those used for Fig. 3 we predict 
a separation line given by: 

(J — K) = 0.34 • (R - K) + 0.19 and R-K > 5.3 (2) 
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Figure 3. Prediction of spectral synthesis models, described in detail in the text. The dashed line is the same in all the panels and 
corresponds to eq. (1). In all the cases, models to the left of the dashed line correspond to ellipticals (thick lines), those to the right to 
dusty star burst (thin lines). 



Not enough observational data are present in this colour- 
colour plane to test the prediction, in particular for dusty 
starburst galaxies; however cluster's galaxies at z=1.26 
(Rosati et al. 1999) and at z=1.27 (Stanford et al. 1997) 
with no emission line fall in elliptical part of R-K vs. J-K 
diagram, while 2 objects belonging to the same clusters but 
with [Oil] emission fall in the starburst region. Indeed at 
this redshift R filter is sampling the rest-frame UV region of 
the spectra and is therefore more sensible than I to residual 
star formation that could affect the colour. I-K colour is less 



affected by this problem, only at higher z (>1.8) ellipticals 
with small residual star formation could fall in the starburst 
region. 

We stress the point that, since z > 2 ellipticals could 
fall in the dusty region (Fig. 1), these colour selection crite- 
ria are valid only in K < 20 — 21 selected sample, where the 
density of z > 2 ellipticals expected in the Pure Luminosity 
Evolution (PLE) models is extremely small (Pozzetti et al. 
1996). Higher redshift ranges can be selected by using dif- 
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Figure 4. Comparison of the selection criterion derived from the 
models with observations. Empty dots are galaxies classified as 
elliptical, solid dots dusty starbursts. When necessary, data were 
transformed into the filter system in use (Cousins' I and Bessel & 
Brett J and K). Data are labelled according to the data source: 
1: Our ERGs selected in HST archive data (Cimatti et al. 1999); 
2: The elliptical in the HDF-South by Stiavelli et al. (1999); 3: 
Ellipticals in a X-ray selected cluster (Rosati et al. 1999); 4: Sub- 
mm sources detected by Smail et al. (1999) (the J magnitude of 
these sources was derived from the published H mag. by assum- 
ing a flat spectrum between these two filters, a very conservative 
assumption); 5 (53W091): Spinrad et al. (1997) and Bruzual & 
Magris (1997); 6 (HR10): Graham & Dey (1996); 7 (HR14): Hu & 
Ridgway (1994); 8: Ellipticals selected by Benitez et al. (1999) in 
th e HDFS with the photometry bv Yahata et al. (200 01 available 
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ferent filters; in particular in the J-K vs. H-K diagram we 
derive the criterion: 

( J - K) = 0.33 • (H - K) + 0.20 and J - K > 2.5 (3) 

to select old (zf > 3) passively evolving elliptical or dusty 
starbursts at 2.0 < z < 2.5. 

In conclusion, the classification of the ERGs using the 
colour-colour diagram I-K vs. J-K seems reliable because 
both model predictions and observed objects give the same 
results. By applying it to large complete samples it will be 
possible to remove most of the foreground stars from the 
ERGs samples and statistically separate ellipticals and star- 
bursts, deriving important information to obtain their for- 
mation and evolution histories. 
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